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Abstract 


We analyze inclusive r decay using a modified version of the a expansion, a 
non-perturbative technique in which the effective coupling is analytic in the 
infrared region. The modification involves renormalization group improve¬ 
ment of the integrand in a spectral representation for the U-function prior 
to implementing the a expansion. The advantage of this approach is that it 
enables us to monitor the structure of the induced power corrections and to 
ensure that these are consistent with the operator product expansion. Nu¬ 
merically the method agrees well with experiment: the comparison is made 
with the physical quantity Rz using Rr as input. 
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I. INTRODUCTION 


Inclusive r~decay provides an ideal arena for extracting information about the low energy 
domain of the strong interactions, in particular the QCD coupling constant as{M^). A 
detailed theoretical analysis of this process has been presented in [|^ (see also There 

has been recent interest in trying to estimate the intrinsic uncertainty in such results due 
to the hnite order truncation of the perturbative series. Attempts at going beyond this 
limitation have included the use of analytic continuation to resum the so-called vr^-corections 


EMU. and infrared renormalon resummations EM (see also PT|-PU|) 
The starting point for the theoretical analysis is the expression E| 


Rt = 


r(r —z/,-hadrons( 7 )) 


r(r- 




'^eh)) 


= 2 



1 + 


( 1 ) 


where R{s) = /mn^y(s -f ie)/7i [] (see e.g. [|^). The integral is not amenable to standard 
perturbation theory as it runs over small values of momentum. However, using Cauchy’s 
theorem, one may rewrite this result as a contour integral in the s-plane with a contour 
running clockwise round the circle |s| = M^, 


R^ = 


1 

2m 


-(1 

z\=l Z 


zf{l + z)D{M^z) 


( 2 ) 


which naively appears to avoid the low momentum region. However, this trick requires the 
Adler D-function, T)(Q^) = —Q‘^dIl{Q‘^)/dQ‘^, to have specihc analytic properties, namely 
to be an analytic function in the Q^-plane except for a cut along the negative real axis. 
This property is broken in perturbation theory due to the presence of the Landau pole at 
= A.QCD- Consequently all perturbative treatments are sensitive to the prescription for 
dealing with the Landau pole. 

m 


For this reason, in a previous publication 


an analysis of r-decay was presented 


making use of an effective running coupling constant which is infrared (IR) hnite and re¬ 
spects the above mentioned analytic properties PT|,|S2[|. While numerically quite successful. 


a conceptual problem with this approach was that one had little knowledge of the structure 
of induced non-perturbative contributions, in particular power corrections. In this letter we 
present a revised formalism which, by hrst analysing the resummmation ambiguity of the 
perturbative series, allows us much greater control over the induced power corrections. In 
Section 2 we explain this modihcation, and in Section 3 apply it to the analysis of r-decay, 
extracting a value of Rz by running the effective coupling up to the Z-scale. Conclusions 
and future applications of this technique are discussed in Section 4. 


II. THE A-EXPANSION FORMALISM 


The nonperturbative expansion method proposed in |3l|j3^ (see also |3^j3^ ) allows one 
to systematically study the low energy regime of QCD and evaluate the integral (j^) either 


^All definitions are as in and we use the standard convention = —q^ where q is the current 
momentum. 
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directly or via use of Cauchy’s theorem as in |]30|j35[| . The method is based on a new 
expansion parameter a connected with ag by the equation 


Air (47r)2 C (1 — a)^ ’ 


( 3 ) 


where C is a positive variational parameter hxed via data from meson spectroscopy |^|. An 


arbitrary Green function may be expanded as a power series in a, a result achieved in practice 
by a resummation of a truncated perturbative series, replacing Ug by a via the appropriate 
expansion of (^). It is clear from (j^) that for all positive values of A the parameter a lies in 
the region 0 < a < 1. 

The evolution of a is dehned by the renormalization group (RG) equation, 


/(a) = /(ao) + 


( 4 ) 


where Oq = cl{QI)- We shall work at O(a^) which, from the expansion of (|^), allows the 


use of perturbative results at O(A^) 
order leads to 


Calculation of the non-perturbative /3-function at this 


/(a) 


1 

5(5+ 3R) 


3 

i=l 


( 5 ) 


where 


J(a, h) = 


2 4 12 

a^b ah‘^ ab 

30 - 216 


9 


(1 - by 


( 1 - 
ln(l — a) 


a)(l — b) 

(2 + by 

63(1 - by 


4 + 126 + 216^ , 

H--In a 


ln(a 


63 

■ 6 ). 


In this expression 


( 6 ) 


where ijk = 123 and cyclic permutations, the values of 6* are the solutions of the equation 
1 + 9a/2 + 2(6 + a)a^ + 5(5 + 3B)a^ = 0, where B = /3i/(2C'/3o), and /3o and /3i are the 
perturbative 1- and 2-loop coefficients of the /3-function. 

For subsequent analysis of r-decay it is useful to separate the higher order terms of the 
D-function via the dehnition 

0(Q")=3SQJ(1 + 4A.„(Q=)), (8) 

t 

where one may consider the contribution to Ag// either from perturbation theory or from the 
a-expansion as appropriate. A direct application of the latter to r-decay, as in [p0| , p5| , then 
corresponds to dehning Ae// via the appropriate expansion of (|^). However, the conceptual 
difficulty with this procedure is that one has little knowledge of the induced non-perturbative 
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corrections contained in the variational series. In order to have control over this aspect we 
must hrst analyse the structure of the RG improved perturbative series. 

To do this in a manner which respects the analytic properties of the D-function we use 
the standard spectral representation, 


r°° 1 

D(Q\X) = Q'‘j 


/o (s + g 2)2 
At next-to-leading-log (NLL) order R{s) is given by |p6 


( 9 ) 


R(s) = 3^Q 


1 + 4A + ( ai - 02 In —1 


( 10 ) 


where 


ai = - [3365 - 22/ - 8C(3)(33 - 2/)], 02 = 4/?o, 


( 11 ) 


It is now convenient to perform an integration by parts in (^, which results in the expression 

1 


roo 

Ae//(t, A)= J druj^r) A +-(oi + 02 - 02 ln(fr))A^ 


( 12 ) 


where t = and we have used (|). The weight function a;(r), given by 


2t 


a; r = 


(l + r)3’ 


(13) 


describes the distribution of virtuality. One may now use RG improvement under the in¬ 
tegral, hrst discussed in |^, with the knowledge that R{s) obeys the same homogeneous 
RG equation as D. Using the two-loop /d-function P{X) = —PoX"^ — PiX^, at NLL order we 
obtain 


roo 

^eff{pX) = d™(r) 


A 


_ ^i^2 ln(l -F/jpAlnr) ^ 


A^ 


l + /doAlnr Po (1-f-/doAlnr)^ (1-1-/dpA Inr)^ 


(14) 


where d = (l/4)(ai -|- 02(1 — In A;)) and A = X{t). This expression is now in a form which 
on inspection will exhibit the divergences associated with large orders in perturbation the¬ 
ory. We observe that the virtuality distribution function (p!3D exactly coincides with the 
function used in ||^ and is numerically very close to that obtained in (see Figure 1), 


which in contrast to the present construction was obtained via an all-orders resummation 
of renormalon contributions in the large Pq limit. 


^InEq. (0), we have restored the factor k representing the renormalization scheme dependence. 
In the so-called V-scheme A:y=l, and in the MS'-scheme which we shall use throughout kjjg = 
exp(-5/3) (see, e.g. @). 
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FIGURES 



FIG. 1. The virtuality distribution functions taken from Ref. (solid line) and the 

function (13) multiplied by a factor of r (dashed line) versus Inr. 


This connection with renormalons is illnstrated more clearly by performing a formal 
Borel transformation on (^), whereby we obtain 


roo 

Ae//(t, A)= / d6exp 
Jo 



B{h), 


(15) 


with 


B{h) = 


1 + 


'A 

2/3o 


+ 2d\h 


r(i + 6/j„)r(2 - hUo) + Ai, 


( 16 ) 


This Borel fnnction exhibits the correct infrared and ultraviolet renormalon poles for the 
H-function, although not the full branch structure [p8| , p9| . Nevertheless, the fact that the 
poles are correctly positioned implies that that the resummation ambiguity associated with 
the hrst IR renormalon is 0{1/Q^)^ which is consistent with the lowest dimension vacuum 
condensate operator in the operator product expansion (OPE) for the H-function. | 

This result is important, as we may now conclude that introduction of the a-expansion at 
this stage, given convergence of the variational series, as has been proven in simpler systems 
EMSB may only induce those non-perturbative power corrections required to cancel this 


^Note that the partial integration performed to obtain dD which generates the appropriate 
virtuality distribution is also responsible for the removal of the first IR singularity at h = l//3o 
required for consistency with the OPE. 

^By convergence we mean strictly absolute convergence of the sequence of approximants for the 
D-function {Di{0{a),Ci), D2{0{a^),C2), ■ ■ ■, T>Ar(0(o^), Gat)}. 
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perturbative ambiguity^. Thus, from the discussion above these corrections will only start 
at 0{1/Q^) and will therefore be consistent with an OPE treatment. 

Introducing the a-expansion for the D-function in this modihed spectral representation 
leads to the replacement of (|T^ by 


fOO _ 

^eff = / dra;(r)A, 
^ 0 


where to O(a^), 


A 


3a^ 

c ^ "c~ ^ c 


d' 


d 

6 + — 

-p- 

10 + 6- 

C 

C 

C 


( 17 ) 


( 18 ) 


The running expansion parameter a = a{Q‘^) is determined via (P, and the variational 
parameter has been determined as O = 21.5 at O(a^) (and also as O' = 4.1 at 0{a^)) 

The H-function constructed in this way may now be utilised in an analysis of r-decay. 


III. r-DECAY ANALYSIS 


Making use of the contour integral representation we isolate the QCD contribution 
by dehning = R^°\l + AR^), with 


= 3(|l4ip + \VJ‘)Sew 


( 19 ) 


The electroweak factor and the CKM matrix elements are Sew = 1-0194, |14d| = 0.9753, 
and 114s I = 0.221 respectively, taken from |l[]. Using the relations (H) and (pTl), and Cauchy’s 
theorem, we obtain^] 


ARr 


48 


rM? ds 

3 Je 




( 20 ) 


which has a modihed kinematic factor as compared to the standard relation (|^) with the 
maximum shifted to s = (2/3)M^. 

Extracting the experimental value of Rr from as described in 0, we obtain = 
3.64 ± 0.02. Using this as input we obtain the effective coupling as{M^) = 0.31 ± 0.01 
at 0{a^). This result is not directly comparable with perturbative extractions of as it 
also includes various non-perturbative corrections, and a correction due to the removal of 
the Landau pole. Nevertheless we note that this quantity is lower than most extractions 
in hxed order perturbation theory, and that this shift is consistent with expectations from 
renormalon resummmations 


®Note that a 1/Q^ correction may also be induced in the running coupling by removal of the 
Landau pole However this has a purely perturbative origin. 

®In this case we take the u,d, and s quarks to be massless and thus ignore threshold effects. 
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It is more consistent to consider only physical qnantities, and therefore we nse the ex¬ 
perimental resnlt for Rt to rnn the conpling np to the Z-scale and extract Rz- In order 
to evalnate this qnantity we apply the matching procednre in the physical timelike channel 
where, at least to leading order, the change in the nnmber of active qnarks is easily associ¬ 
ated with the threshhold for pair prodnction. The effective conpling and its derivative are 
reqnired to be continnons at the threshold points, which is acheived by solving the continn- 
ity eqnations for Cf and a^. Using the standard heavy qnark masses, rric = 1.6 GeV and 
rrib = 4.5 GeV, and applying this matching procednre we hnd 

= 20.96 ± 0.01 at 0{a^), (21) 

= 20.83 ±0.01 at O(a^), (22) 


which, at 0{a^), is within one standard deviation of the experimental resnlt Rz = 20.77±0.07 
Note that this is a considerably better £t to the data than using the technique of 


30 


when one accounts for the change in the data over the intervening period. 


IV. CONCLUDING REMARKS 


In the present letter we have described a development of the technique introduced 
in which allows significantly greater control over the stucture of induced non- 

perturbative corrections. This has allowed an application to r-decay and the resulting 
normalisation of the effective coupling results in very good agreement with experimental 
data at the Z-scale. 

The intriguing possibility raised by this technique is that, as discussed in section 2, 
provided one assumes convergence of the variational series, the induced power corrections 
will be consistent with the existence of vacuum condensate operators in the OPE. Thus one 
may wonder whether this formalism may allow an investigation of the meson spectrum in 
the framework of QGD sum rules which does not explicitly include arbitrary condensate 
parameters (see e.g. [^,|^), the assumption being that these corrections, at least in some 
averaged sense, are automatically induced by convergence of the series. The results of an 
investigation of this kind will be presented in a forthcoming publication . 


The authors thank D.V. Shirkov for discussions, and the support of A.R. by the Gom- 
monwealth Scholarship Gommission and the British Gouncil, and the partial support of 
I.L.S. by a Royal Society Grant is gratefully acknowledged. 
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